Fermions, strings, and gauge fields in lattice spin models 
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We investigate the general properties of lattice spin models with emerging fermionic excitations. 
We argue that fermions always come in pairs and their creation operator always has a string-like 
structure with the newly created particles appearing at the endpoints of the string. The physical 
implication of this structure is that the fermions always couple to a nontrivial gauge field. We 
present exactly soluble examples of this phenomenon in 2 and 3 dimensions. Our analysis is based 
on an algebraic formula that relates the statistics of a lattice particle to the properties of its hopping 
operators. This approach has the advantage that it works in any number of dimensions - unlike the 
fiux-binding picture developed in FQH theory. 
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I. INTRODUCTION 

For many years, it was thought that Fermi statistics 
were fundamental, in the sense that one could only ob- 
tain a theory with fermionic excitations by introducing 
them by hand (via anti-commuting fields). Then, over 
the past two decades, this view began to change. A 
number of real world and theoretical examples showed 
that fermions and anyons could emerge as low energy 
collective modes of purely bosonic systems. The first 
examples along these lines were the fractional quantum 
Hall states. [1, 2] Usually we think of the FQH states as 
examples of anyonic excitations emerging from interact- 
ing fermions. However, from a purely theoretical point of 
view, the same effects should occur in systems of interact- 
ing bosons in a magnetic field. [3] There were also indica- 
tions of emerging fermions in the slave-boson approach 
to spin- 1/2 systems [4-10] and in the study of resonat- 
ing valence bond (RVB) states. [11-14] Unfortunately, the 
RVB picture and the slave boson approach both rely on 
approximate or mean-field techniques to construct and 
analyze these exotic states. More recently, a number 
of researchers have introduced exactly soluble or quasi- 
exactly soluble models with emerging fermions. [14-21] 
These models allow for a more well-controlled analysis, 
albeit in very specific cases. 

The mean-field approach and the exactly-soluble ex- 
amples both provide clues to the structure and basic 
properties of bosonic models with fermionic excitations. 
They indicate, among other things, that fermions never 
appear alone in lattice spin systems. Instead they al- 
ways come together with a nontrivial gauge field and 
the emerging fermions are associated with the deconfined 
phase of the gauge field. [7, 9, 10] The deconfined phases 
always seem to contain a new kind of order - topological 
order, [10, 22] and the emerging fermions and anyons are 
intimately related to the new order. 

This was particularly apparent in the context of the 
slave-boson approach. In this technique, one expresses 
a spin- 1/2 Hamiltonian in terms of fermion fields and 
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gauge fields. [4, 5] Clearly, the presence of fermion fields 
does not, by itself, imply the existence of fermionic quasi- 
particles - the fermions appear only when the gauge field 
is in the deconfined phase. Ref. [6-10] constructed several 
deconfined phases where the fermion fields do describe 
well defined quasi-particles. Depending on the properties 
of the deconfined phases, these quasi-particles can carry 
fractional statistics (for the chiral spin states) [6-8] or 
Fermi statistics (for the Z2 deconfined states). [9, 10] 

Although it was less evident initially, a similar picture 
has emerged from the study of RVB states and the as- 
sociated quantum dimer models. Ref. [11, 12] originally 
proposed that fermions (spinous) could emerge from a 
nearest neighbor dimer model on a square lattice. Later, 
it was realized that the fermions couple to a Z2 gauge 
field, and the fermionic excitations only appear when the 
field is in the deconfined phase. [13] It turns out that the 
dimer liquid on square lattice with only nearest neighbor 
dimers is not in the deconfined phase (except at a crit- 
ical point). [12, 13] However, on a triangular lattice, the 
dimer liquid does have a Z2 deconfined phase. [14] The re- 
sults in Ref. [11, 12] are valid in this case and fermionic 
quasi-particles do emerge in a dimer liquid on a triangu- 
lar lattice. 

In this paper, we attempt to clarify these observations 
and to put them on firmer foundations. We give a gen- 
eral argument which shows that emerging fermions al- 
ways occur together with a nontrivial deconfined gauge 
field. In addition, we derive an algebraic formula that 
allows one to calculate the statistics of a lattice particle 
from the properties of its hopping operators. We feel that 
this formula both elucidates the fundamental meaning of 
statistics in a lattice system and simplifies their compu- 
tation. This approach has the added advantage that it 
works in any number of dimensions - unlike the flux bind- 
ing picture [3] developed in FQH theory. We would like 
to point out that all the previously mentioned examples 
of emerging fermions are two dimensional models. The 
emerging fermions in those models are related to the flux 
binding picture in Ref. [3] . Our algebraic approach allows 
us to establish the emergence of fermions in an exactly 
soluble 3D bosonic model. 

Our paper is organized as follows: in the first section, 
we give a definition of the statistics of a lattice particle. 
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In the second section, we derive the algebraic statistics 
formula discussed above. In the third section, we apply 
the formula to the case of fermionic (or anyonic) excita- 
tions in a lattice spin system. The formula demands that 
the fermions always come in pairs, and that their pair cre- 
ation operator has a string-like structure with the newly 
created particles appearing at the ends. We show that 
the strings represent gauge fluctuations, and the fermions 
carry the corresponding gauge charge. After presenting 
the general argument, we devote the last two sections to 
exactly soluble examples of this phenomenon in 2 and 3 
dimensions. 

The above result suggests a very interesting picture of 
emerging fermions. Fermions and gauge fields appear to 
be two sides of the same coin: in some sense, fermions are 
the ends of strings and gauge fields are the fluctuations 
of strings. Extended string- like structures seem to be 
the key to understanding both Fermi statistics and gauge 
fields. 



II. DEFINITION OF STATISTICS 

What do we mean when we say that a particle is a 
"boson" or a "fermion" ? The most common way to define 
statistics is to use the algebra of creation and annihilation 
operators. If these operators satisfy the bosonic algebra, 
the corresponding particles are said to be bosons; if they 
satisfy the fermionic algebra, the particles are said to be 
fermions. 

In this paper, we will use another, equivalent, defi- 
nition based on the path integral formulation of quan- 
tum mechanics. Given any m particle system with short- 
range interactions, we can write down the corresponding 
multi-particle action S. The most general action S has 
two components: 

S = Sloe "I" Stop 

The first term, Sioc, can be any local expression. Typ- 
ically, it's just the usual classical action. The second 
term, Stop, is less familiar. In general, it can be any 
expression which only depends on the topology (e.g. ho- 
motopy class) of paths. In 3 or more dimensions, the 
form of Stop is highly constrained. The amplitude for a 
closed path P must be of the form 

where n is the number of particle exchanges that occur 
in P. The statistics of the particles are defined by the 
sign in this expression. If this sign is positive, we say 
the particles are "bosons", if the sign is negative, the par- 
ticles are "fermions". In 2 dimensions, there are other, 
more complicated possibilities for Stop- These particles 
are called "anyons." 

According to this definition, the statistics are com- 
pletely determined by the topological term Stop in the 
action. One way to isolate this term is to compare the 
amplitude for two paths which are the same locally (this 
will be made more precise in the next section) but differ 



in their global properties. In particular, we can com- 
pare the amplitude for a path which exchanges two par- 
ticles with another path which is the same locally, but 
doesn't exchange the two particles. The difference be- 
tween the phases of the two amplitudes is then precisely 
giStop _ (-1-1)^ the statistical phase. 

We now reformulate this (theoretical) test of statistics 
so that it can be applied to particles on a lattice. The 
prescription is as follows: we take a two particle state 
|ri,si), and consider a product of hopping amplitudes 
along a lattice path which exchanges the particles 

(r„, s„|ff |r„_i, s„_i)(r„_i, s„_i|if |r„_2, Sn-2) 

■ ■■{r3,S3\H\r2,S2){r'2,S2\H\ri,si) (1) 

Each hopping amplitude in this expression involves one 
of the two particles moving to a neighboring site while 
the other particle remains fixed. At the end, the two 
particles have exchanged places: r„ — si, s„ = ri. 

Now we construct another path which is the same lo- 
cally, but doesn't exchange the two particles, and take 
a product of hopping amplitudes along this path. The 
claim is that the difference between the phases of these 
two expressions is precisely the statistical phase of the 
particles. That is, the particles are bosons, fermions, or 
anyons depending on whether the phase difference is -|-1, 
— I, or something else. 

One way to see this is to remember the derivation of 
the path integral formulation of quantum mechanics. Ac- 
cording to the standard derivation, the amplitude e'"^ for 
a two particle path \r{t), s{t)) is given by a product 

e'^ oc (r(i„),s(i„)|e-^^^*|r(t„_i),s(t„_i)) 
...{r{h),s{h)\e-^"^*\r{t2),s{t2)) 
(r(i2),s(i2)|e-^^^*|r(ti),s(ii)) 

in the limit that At = t„ — t„_i 0. 

Now, in the discrete (lattice) case, the above expression 
can be further simplified. Since At — > 0, we can rewrite 
it as 

e'^ (X {r{tn),s{tn)\{l - iHAt)\r{tn-i),s{tn-i)) 
...{r{t3),s{hm-iHAt)\r{t2),s{t2)) 
{rit2),s{t2m-iHAt)\r{h),s{h)) 

Successive states |r(fi), s(ti)), |r(ti_|_i), s(ti_|_i)) must ci- 
ther be identical, or must differ by a single particle hop. 
The matrix elements between identical states don't con- 
tribute to the phase. Thus, we can drop them without 
affecting our result. We are left with: 

e'' oc {r{t',),s{t',)\{-tHAtMt',_,),s{t',_,)) 
...{rit',),s{m-zHAt)\rit',),sit',)) 
{rit',),s{t',)\i-zHAt)\rit[),s{t[)) 

where the t^s are all distinct. 

Now, as we discussed earlier, the statistical phase can 
be obtained by comparing the phase of this amplitude 
with another amplitude which is the same locally but 
doesn't exchange the two particles. When we make this 
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comparison, the factors of — i drop out (since they con- 
tribute equally to the two products). Thus, it suffices to 
compare products of the form: 

(r(4),s(4)|i7|r(4_i),s(4_i)) 
...{rit',),sit',)\H\r{f',),.s(f',)) 
{r{t',),s{t',)\H\rit[),s{t[)) 

This is precisely the expression in Eq. (1). 

III. STATISTICS AND THE HOPPING 
OPERATOR ALGEBRA 

In this section we derive a simple algebraic formula 
for the statistics of a particle hopping on a lattice. This 
formula is completely general and holds irrespective of 
whether the particles are fundamental or are low energy 
excitations of an underlying condensed matter system 
(e.g. quasi-particlcs) . 

We begin with a Hilbert space which describes n hard- 
core particles hopping on a d dimensional lattice. The 
states can be labeled by listing the positions of the n 
particles: \ii,i2, ■■■in)- The particles are identical so 
the states \ii,i2, ■■■in) do not depend on the order of 
ii,i2, ■■■,in- For example, 

\ii,i2, ■■■in) = \i2,ii, ■■■in) 
A typical Hamiltonian for this system is of the form 

<ij> 

where tij are "hopping operators" with the property that 

tij\j,ii,^^;in-i) oc ...,i„_i) (2) 

We assume that the hopping is local, i.e 

[tij,tki]=0 (3) 

if i,j,k,l are all different. 

Our goal is to compute the statistics of the particles 
described by this hopping Hamiltonian H. In the fol- 
lowing, we show that the statistical angle can be derived 
from the simple algebraic properties of the hopping op- 
erators. More precisely, we show that the particles obey 
statistics e'^ if 

til^kitij — ^ iijiki^il (4) 

for any three hopping operators tij,tki,tii, where j,k,l 
are (distinct) neighbors of i (ordered in the clockwise 
direction in the case of 2 dimensions). The orientation 
convention in 2 dimensions is necessary for the anyonic 
case. 

The simplest examples where we can apply this for- 
mula are the cases of noninteracting hard-core bosons or 
fermions. In these cases the n particle Hamiltonian can 

be written as 

<ij> 



(a) (b) 




j 5 j' j 5 j' 

FIG. 1: (a) The path of the particles in the construction of the 
swapped state swapped). The numbers label the order in 
which the paths are traversed, (b) The path of the particles 
in the unswappcd state unswapped). Notice that the two 
states only differ in the ordering of the paths. However, in 
one case the particles switch places, while in the other, they 
do not. 

where the Cj's are the boson or fermion annihilation op- 
erators. The hopping operators are just Uj = —tcjcj. A 
little algebra confirms that the boson and fermion hop- 
ping operators do indeed satisfy Eq. (4) with e'^ = 
and —1, respectively. 

We now give a general derivation of the formula. We 
begin with the state |i, ...) which contains particles at 
sites i,j and other particles far away. 

Imagine that we exchange the two particles at i,j us- 
ing an appropriate product of hopping operators. One 
expression for the final swapped state is 

swapped) = {tjji){tj>p...tqr){tiltim--tnj') 
{tj'j){trs^^^tti)\i,j,^^^) 

Here, we've grouped the hopping operators into 5 terms. 
Each of these terms moves a particle along some path 
on the lattice. The cumulative effect of all 5 terms is to 
exchange the particles (see Fig. la). 

Similarly, we can construct a path which doesn't swap 
the two particles. One expression for the final unswapped 
state is 

unswapped) = {tjj>){tj,j){tiitim...tnj') 
(ij/p...ig7.) (^7-5 ...i^i) j^, J, ...) 

Once again, we've grouped the hopping operators into 5 
terms. In this case, the cumulative effect of these terms 
is simply to move the two particles along independent 
loops without any swapping taking place (see Fig. lb). 

Intuitively, we expect that the phase difference be- 
tween the swapped and unswapped states is related to the 
statistical phase of the particles. In fact, the special way 
we've constructed our states allows us to make a much 
stronger statement. Notice that the two states involve 
the same product of hopping operators. The only differ- 
ence is the order of these operators. This means that, 
in the swapped and unswapped states, the two parti- 
cles trace out the same total path in configuration space. 
This is important because it implies that any phase which 
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comes from gauge fields or other Berry phases contributes 
equally to the swapped and unswapped states. The phase 
difference between the swapped and unswapped states is 
therefore exactly equal to the statistical phase of the par- 
ticles. 

This intuitive argument can be made rigorous using the 
discussion in the previous section. Indeed, it's not hard 
to see that the two space-time paths traced out by the 
swapped and unswapped states are exactly the same lo- 
cally - they only differ by a reordering. Furthermore, the 
phase difference between the swapped and unswapped 
states can be written as a phase difference between two 
expressions of the form (1). Thus, taking the phase dif- 
ference between the swapped and unswapped states is 
completely equivalent to the procedure derived in the 
previous section. 

We now compute this phase difference. We use the 
assumed algebraic relation 

ij'pinj'ij' j ~ ^ ^j' j^Tij'^j'p 

Applying this relation to the swapped state and reorder- 
ing the hopping operators using locality, we find 

swapped) = e'^ |i, j, unswapped) 

The phase difference is e*^, so the particles obey statis- 
tics e*^, as claimed. A similar formula can be derived for 
the relative statistics of distinguishable particles. This 
formula, together with its derivation, is given in the ap- 
pendix. (See Eq. (Al)). 

IV. FERMIONS AND STRINGS 

In this section we consider the properties of a lattice 
spin system with fermionic or anyonic excitations. We 
argue that the excitations are always created in pairs, and 
the creation operator for a pair of particles has a string- 
like structure, with the new particles located at the ends. 
One interpretation of this is that fermions never appear 
alone - they always come with some kind of gauge field. 

We now make these statements more precise. Suppose 
we have a lattice spin system with fermionic excitations 
(the anyonic case is completely analogous). The total 
Hilbert space of the lattice spin model is a direct prod- 
uct of local Hilbert spaces Hi, each associated with an 
individual lattice site /. We expect that the total Hilbert 
space contains a low energy subspacc which is spanned by 
n fermion states. The corresponding low energy effective 
Hamiltonian is given by 

iJeff = PHP 

where P is a projection operator onto the n fermion sub- 
space. Typically, iJeff can be written as a sum 

HeS = ^ {tij + tji) 
<ij> 

where the Uj are hopping operators. We expect that the 
tij are local in the underlying spin degrees of freedom. 



That is, they act only within the local Hilbert spaces 
near sites i, j. [30] We would like to point out that the 
lattice on which the fermion hops (labeled i) may not be 
the same as the lattice formed by the local Hilbert spaces 
(labeled /). 

Now, consider the string-like product of hopping oper- 
ators 

Wir — tijtjfztfzl...tpqtqy 

This operator destroys a fermion at site r and creates one 
at site i: 

Wir\r,...) oc 

where |r, ...), are low-energy states with fermions 

at r, i, (and possibly other fermions far away). Further- 
more, this operator clearly has a string-like structure: it 
is made up of a string-like product of operators, each 
of each is local in the underlying spin degrees of free- 
dom. The only issue is that this string operator might 
be trivial. That is, Wir might act trivially on all the in- 
termediate spins at j, k, q, and only have a nontrivial 
effect near i and r. In that case, Wir isn't really a string 
at all. 

This is a legitimate concern, since the string operator 
typically is trivial in the case of bosonic quasi-particles. 
Consider, for example, a lattice of noninteracting spins 
in a magnetic field: H = —BJ2i<^i- The ground state 
has af = 1 for all i. The excitations are obtained by 
flipping one spin: (t| = —1 for some j. If we perturb the 
Hamiltonian by a term t^^^j^ alaj , these excitations 
acquire dynamics. The corresponding hopping operators 
are then tij = ^crluj and the string operator is 

^Vir — tijtjk-'-tqr 

OC (aiaj)(a]ai)...(a>i) 
= ajal. 

The operator Wir clearly creates a particle at position i 
and destroys one at position r. However, as we see above, 
it is trivial and has no string-like structure. The particle 
creation and annihilation operators are completely local. 

Our claim is that this can never happen in the case of 
emerging fermions. That is, the string operator Wir can 
never be written as a product 

Wir = AiBr 

where Ai, Br are operators which act only within the 
local Hilbert spaces near sites r. Physically, this means 
that the fermion creation and annihilation operators are 
never local in the underlying bosonic degrees of freedom. 
Instead, they naturally appear in pairs, with a string-like 
structure connecting them. [31] 

Proof: We wish to use the algebraic formula (4) from 
the previous section. That formula relied on the locality 
of the hopping operators (3). Here we would like to use 
a weaker locality condition: 

[tij,tki\=^ (5) 



if i, j, are far from k, I. One can make a similar argument 
using this weaker assumption. One arrives at a slightly 
weaker algebraic relation: 

WuWkiWij = -WijWkiWa 

if i, I, k are sufBciently far from each other. Since 
WkiWij = Wkj, we can rewrite this as 

WuWkj = -WijWki (6) 

Condition (6) is a direct consequence of the fermionic 
nature of the quasi-particles. We will now show that it 
implies that the strings are nontrivial. Indeed, if the Wir 
could be written as Wir = AiBr, then 

WiiWkj = AiBiAkBj 
WijWki = AiBjAkBi 

But the v4's and _B's arc local operators, so they commute 
with each other when they are well separated. We can 
therefore rearrange the operators in these two equations 
to obtain 

WuWkj = WijWki 

This directly contradicts the fermion condition (6). We 
conclude that the string operator is always nontrivial. □ 
The presence of this nontrivial string indicates that 
fermions always appear together with some kind of gauge 
field. One way to see this is to consider a closed loop 
of hopping operators tijtjk---tpqtqi. This string can be 
interpreted as a Wilson loop operator, since its phase is 
precisely the accumulated phase of the particle when it 
traverses a loop. The fact that it is nontrivial (that is, not 
equal to the identity operator) means that the particle is 
coupled to a nontrivial gauge field. 



V. A 2D EXAMPLE 

In this section, we present an exactly soluble lattice 
spin model with fermionic excitations. The exactly sol- 
uble model provides a concrete realization of the string 
picture discussed above. The emerging fermions turn out 
to be coupled to a Z2 gauge field. 

In this spin-1/2 system, proposed by Kitaev[15], the 
spins live on the links of a square lattice (Fig. 2). The 
Hamiltonian is 



(7) 



Here i labels the links. I labels the sites, and p labels 
the plaquettes of the square lattice. Also, Cj denotes the 
loop connecting the four spins adjacent to site I, while 
Cp denotes the loop connecting the four spins adjacent 
to plaquette p (see Fig. 2). 

This model is exactly soluble since all the terms in the 
Hamiltonian commute with each other. The ground state 
satisfies 




T , t , t T 



FIG. 2: A schematic diagram depicting Kitaev's model. Each 
spin is drawn as a dot. The paths C'j and Cp arc drawn as 
dotted diamonds, and are labeled with I's or 3's according to 
whether the corresponding term in the Hamiltonian involves 



for all sites /, and 



for all plaquettes p. There are two types of (localized) 
excited states. We can have a site where 



n 



a\ = -1 



or we can have a plaquette where 



-1 



We call the first type of excitation a "charge" and the 
second type of excitation a "flux." Static charge and flux 
configurations are exact eigenstates of the above Hamil- 
tonian. Thus, the charge and flux quasi-particles have no 
dynamics. 

This lack of dynamics is a special feature of the above 
model. However, we arc interested in the properties of a 
generic Hamiltonian in the same quantum phase. Thus, 
we need to perturb the system, and analyze the resulting 
dynamics. The simplest nontrivial perturbation is 



(8) 



It's not hard to see that the first term allows the fluxes 
to hop from plaquettes to adjacent plaquettes, while the 
second term allows the charges to hop from sites to neigh- 
boring sites. 

We now calculate the statistics of the fluxes. To do 

this, we restrict our Hamiltonian to the low energy sub- 
space with n fluxes and zero charges. Within this sub- 
space, our Hamiltonian reduces to 



To make contact with our previous formalism, we write 
this as 

<pi> 
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where the sum is taken over adjacent plaquettes, p, q, tpq 
is defined by tpq = ^cr|, and i the hnk joining p and q. 

To calculate the statistics, we need to compare 
tpqtrptps with tpstrptpq. WcU, it's obvious from the def- 
inition that all the hopping operators tpq commute with 
one another. Therefore, 



tpqtjtptpg 



tpstrptpq 



We conclude that the fluxes are bosons. In the same way, 
one can show that the charges arc also bosons. 

Next, we consider the bound state of a flux and a 
charge. That is, we consider excitations with a flux 
through some plaqucttc p and also a charge at one of 
the sites / adjacent to p. 

These bound states are not actually stable for the 
above Hamiltonian (8) - the charge and flux will separate 
from one another over time. However, one can imagine 
modifying the Hamiltonian so that charges and fluxes 
prefer to be adjacent to each other. In this case, the 
bound state is a true quasi-particle. 

Suppose we've made such a modification. We can then 
consider the statistics of the bound state. If we restrict 
our Hamiltonian (8) to the low energy subspace with n 
bound states, then our Hamiltonian reduces to 



To understand the effect of these terms, imagine we have 
a bound state with a flux at p and a charge at I. It's not 
hard to see that the flrst term allows the flux to hop to the 
two neighboring plaquettes which are also adjacent to /. 
Similarly, the second term allows the charge to hop to the 
two neighboring sites, adjacent to p. All other hopping 
destroys the bound state and is therefore forbidden for 
energetic reasons. Thus, the perturbation gives rise to 
4 types of hopping operators - 2 corresponding to fluxes 
and 2 corresponding to charges. (See Fig. 3). 

Formally, we can write our low-energy Hamiltonian as 



X] {i{p,I)(q,J)+i{q,J){p,I)) 
<(p,I)(q,J)> 



where the hopping operators are defined by 



t{p,I){q,J) 



■hi J3 3 
y^i or -a, 



depending on whether (p, /), (q, J) differ by a flux hop 
or a charge hop. In the flrst case, i is deflned to be the 
link joining p and q, while in the second case, i is the 
link joining / and J. 

To calculate the statistics we need to compare a prod- 
uct of the form tit2t3 with the product t;3,t2ti, where the 
tk are hopping operators involving a single bound state 
at position (p, /) . (We could write out these expressions 
precisely, but it's messy and not very enlightening). 

Now, as we discussed above, there are four different 
ways that a bound state at (p, I) can hop - two charge 
hops, and two flux hops. Each of the flux hopping op- 
erators can be paired with a corresponding charge op- 
erator which involves the same link i. It's easy to see 




FIG. 3: A bound state of a charge, denoted by a dark circle, 
and a flux, denoted by a shaded square. The arrows illustrate 
the four ways the bound state can hop. The two charge hops 
are denoted by solid arrows, and have hopping operator a^. 
The flux hops are denoted by dotted arrows, and have hopping 
operator . Notice that each charge hopping operator anti- 
commutes with a corresponding flux hopping operator, but 
everything else commutes. 



that each of the flux operators anti-commutes with the 
corresponding c:liarg(^ operator (since one involves a cr^, 
while the other involves a cr"^). However, everything else 
commutes. (See Fig. 3). 

With these facts in mind, we can compare tit2i3 with 
tzt^ti- By our discussion above, any set of 3 hopping 
operators involving (p, I) must contain exactly 2 which 
anti-commute. Thus, exactly 2 of ti,t2,t3 anti-commute. 
This implies that 



tlt2t3 = —t3t2tl 

We conclude that the bound states are fermions. Of 
course, this result is not that surprising once we no- 
tice that the charges and fluxes have relative statistics 

giljirel. — ^ 

We can easily check that the charges and fluxes have 
relative statistics —1. Normally, we would just use the 
formula (Al) derived in the appendix. However, this 
relation assumes that the two types of particles both hop 
on the same lattice. In our case, the charges hop on the 
square lattice, while the fluxes hop on the dual lattice. 

It's not hard to see that in our case the formula (Al) 
needs to be modified to 



tijtf 



ipq^IJ 



where J, J are neighboring sites, p, q are neighboring 
plaquettes, and the links connecting I, J, and p, q are 
the same. From our previous calculations we know that 
tpq = ^crl, and ijj = ^cr?. Therefore, tpq, tjj anti- 
commute, and the relative statistics are e"'^'"' = — 1. 

We now see why the two particles were called "charges" 
and "fltixes" - they have the same statistics and relative 
statistics as Z2 charges and fluxes. It turns out that this 
connection with Z2 gauge theory extends beyond the low 
energy regime - in fact, all the way to the lattice scale. 
One can show that the Kitaev model is exactly equivalent 
to standard Z2 gauge theory coupled to a Z2 Higgs field. 

We argued earlier that whenever fermions or anyons 
occur in a bosonic system, they are always created in 
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FIG. 4: Three examples of string-like creation operators. The 
string operators are drawn together with the particles they 
create at their endpoints. The thick solid line (drawn on 
the lattice) is an example of a charge string operator. The 
dotted line (drawn on the dual lattice) is an example of a 
flux string operator. At the bottom, we give an example of a 
string operator for the bound state of a charge and flux. It is 
essentially a combination of a charge string and a flux string. 



the charge and flux strings. Let P = (pj^, 7i)...(p„, 7„) 
be a path in bound state configuration space. Then the 
associated bound state string operator is 

W{P)'xYlat: (11) 

k 

where are the links along the path and 0^ = 1 

or 3 depending on whether {pi..,Ik) differ by a flux hop 
or a charge hop, respectively (see Fig. 4). Once again, 
one can show that the string W{P) is a creation operator 
for a pair of bound states located at the ends of P. 

In each of these examples, the creation operator for 
fermions or anyons has an extended string-like structure. 
It is important to note that the position of this string 
is completely unobservable. That is, in each case, the 
excited state VF(P)|0) is independent of the position of 
P: it only depends on the position of the endpoints of P. 



VI. EXACTLY SOLUBLE 3D MODEL 



pairs, and the pair creation operator has a string- like 
structure. The above exactly soluble model (7) provides 
a good example of this phenomenon. 

We begin with the charges. We can construct the 

string operators associated with these particles by tak- 
ing products of their hopping operators along some path 
P = Ii, ...In on the lattice. We find 



W{P) 



oc 



.tl.I..-.tT 



(9) 



where i\...in-i are the links along the path (see Fig. 4). 
Notice that W{P) is nontrivial - it acts non- trivially on 
the spin degrees of freedom along the string P. 

It is also a pair creation operator: If we apply W{P) 
to the ground state, the resulting state is an exact eigen- 
state with two charges - one located at each endpoint of 
P. One way to see this is to notice that W{P) commutes 
with everything in the Hamiltonian except and 

Wc'j uj. The string anti-commutes with these two op- 
erators, so when we apply it to the ground state, we get 
a state with two charges located at Ii,In- 

The case of the flux quasi-particlcs is very similar. We 
take the product of the flux hopping operators along some 
path P' = Pi, ...Pn on the dual lattice. We find 



W'{P') 



''Pn-2Pn.-l^Pn- 



(10) 



where ii...z„_i are the links along the path (see Fig. 4). 
Just as before, one can show that the string W'{P') is 
a creation operator which creates two fluxes at the end- 
points of P' . 

Finally, consider the case of the bound state of the 
charge and the flux. The hopping operator for the bound 
state is a combination of the charge and flux hopping 
operators, so the string turns out to be a combination of 



In two dimensional systems, one can create a fermion 
by binding a Z2 vortex to a Z2 charge. This is how we 
obtained the fermion in the above spin- 1/2 model on the 
square lattice. Since both the Z2 vortex and Z2 charge 
appear as the ends of open strings, the fermions also ap- 
pear as the ends of strings. However, in three dimensions, 
we cannot change a boson into a fermion by attaching a 
TT-flux. Thus one may wonder if fermions still appear as 
the ends of strings in (3-|-l)D. In this section, we study an 
exactly soluble spin-3/2 model on a cubic lattice. We will 
show that the creation operator for fermions does indeed 
have a string-like structure. This example demonstrates 
that the string picture for fermions is more general then 
the flux-charge picture. 

Our model has four states for each site of a cubic 
lattice. Thus we call it a spin-3/2 model. Let 7"**, 
a, 6 e {x, X, y, y, z, z}, a 7^ &, be 4 x 4 hermitian matrices 
that satisfy 

j«-b — _jba_ (^a6-jt 

[y"\ 7'=''] =0, if a, b, c, d are all different 
7«Y^=i7«^, ajic 



(12) 



A solution of the above algebra can be constructed by 
taking pairwise products of Majorana fermion operators 

A", a e {x,x,y,y,z,z}: 



{A",A''}=2<5„6 



(13) 



The six Majorana fermion operators naturally require 
a space of dimension 2^/^ = 8, but if we restrict them to 
the space Ha A" = 1, we obtain the desired 4x4 hermi- 
tian matrices. Alternatively, a more concrete description 
of the 7"'' is given the appendix, where we express the 
7"'' in terms of Dirac matrices. 



In terms of 7?'', the exactly soluble spin-3/2 Hamilto- 
nian can be written as 

H^^gJ^Fp (14) 

P 

where p labels all the square plaquettes in the cubic lat- 
tice, and the Fp arc "flux" operators defined by 

rp yx xy yx xy 

zy yz zy yz 
or ^'y . ^ ^'^f 

Ji h+y h+y+z H+z^ 

or iriU^tu+^Z^- (15) 

depending on the orientation of the plaquette p. Just 
as in the Kitaev model, [15] all the Fp commute with 
each other and all the Fp have only two eigenvalues: ±1. 
Thus, we can solve the Hamiltonian by simultaneously 
diagonahzing all the Fp. If |{/p}) is a common eigen- 
state of Fp with Fp\{fp}) = /p|{/p}), fp = ±1, then it 
is also an energy eigenstate with energy 

E{{fr>}) = -gY.fv 
p 

The one subtlety is that the fp are not all independent. If 
C is the surface of a unit cube, then we have the operator 
identity 

n^p-1 

pec 

It follows that 

n/p=i 

pec 

for all cubes C. This constraint means that the spectrum 
of our model is identical to a Z2 gauge theory on a cubic 
lattice. The ground state can be thought of as a state 
with no flux: fp — 1 for all p. Similarly, the elementary 
excitations are small flux loops where 

/p--l 

for the four plaquettes p adjacent to some link {ij). We 
can think of these excitations as quasi-particles which live 
on the links of the cubic lattice. 

We would like to compute the statistics of these exci- 
tations. It is tempting to assume that they are bosons 
since the model is almost the same as a Z2 gauge the- 
ory. However, as we will show, this superficial similarity 
is misleading: the flux loops are actually fermions. 

As in the Kitaev model, the quasi-particles are exact 
eigenstates of our Hamiltonian. Thus, they have no dy- 
namics and it is difficult to compute their statistics. How- 
ever, this is a special feature of our model. We need to 
perturb the theory to understand generic states in the 
same quantum phase. The simplest perturbation is 
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FIG. 5: A fermion on link 1 can hop onto 10 different links; 
2-11. The hopping 1 — > 2 is generated by 7?^^, the hopping 
1 10 by 7?^, and the hopping 1 ^ 4 by 7^^. 

To compute the statistics, we need to restrict ourselves 
to the low energy subspace with n quasi-particles. In this 
subspace, H' reduces to 

i.a.h 

The effect of each term 7?^ is to allow the quasi-particles 
to hop between the two links -\- a)) and + b)) 
adjacent to site i. Here a — x ii a = x, a = —x if a = x, 
etc . (See Fig. 5.) Thus, our Hamiltonian can be written 
in the standard hopping form: 

■f^off = Y{t(i[i+a.)){i{i+b)) +^{i(i+b)){i(i+a))) 
i,a,b 

where t(i(i+a)>(i(i+6)> = bt'- 

To calculate the statistics we need to compare a prod- 
uct of the form tit2t^ with the product t^t2ti, where the 
tk are hopping operators involving a single link {ij). (As 
in the Kitaev model, we could write out these expressions 
explicitly, but it's not very enlightening). 

Note that a quasi-particle on the link {ij) can hop to 
any of the 10 neighboring links, 5 of which are adjacent 
to site i, and 5 of which are adjacent to site j (see Fig. 5). 
Using the algebra (12) or the Majorana fermion represen- 
tation, we find that the 5 hopping operators associated 
with site i all anti-commute with each other and simi- 
larly for site j . On the other hand, each of the operators 
associated with i commute with each of the operators 
associated with j . 

With these facts in mind, wc can now compare tit2i3 
with ^3^2^!- There are essentially two cases: either all 
three of the i^'s are associated with a single site i or j, 
or two involve one site and one involves the other. In 
the first case, all the tk's anti-commute; in the second 
case, a pair of the t^'s anti-commute and everything else 
commutes. In either case, we have 

We conclude that the quasi-particles are indeed fermions. 
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FIG. 6: An example of a string in the 3D model. 
This string corresponds to the operator W{P) oc 



(a) 



(b) 



5; 



FIG. 7: (a) The path of the particles in the construction of 
the wound state wound). The dotted line is the path of 
particle 1, and the solid line is the path of particle 2. The 
numbers label the order in which the paths are traversed, (b) 
The path of the particles in the unwound state unwound) . 
Notice that the two states only differ in the ordering of the 
paths. However, in one case the particles wind around each 
other, while in the other, they do not. 



Next, we construct the associated string operator. 
Taking a product of hopping operators along a path (of 
Hnks), P = (^1*2), (^2*3), ■■■{in~iin), gives 

WiP)^jt:''^ltf^...7t" (16) 

where &„j, a„j+i is the pair of indices associated with the 
link {imim+i}- (See Fig. 6.) Notice that the string opera- 
tor commutes with the Hamiltonian, except at the ends. 
Thus, when we apply W{P) to the ground state, the only 
effect is to flip the signs of fp for the plaquettes p adja- 
cent to the links {1112), {in-iin)- This means that the 
open string operator creates a pair of fermions at its two 
ends. 



VII. CONCLUSION 

In this paper we have derived a general relation (4) be- 
tween the statistics of a lattice particle and the algebra of 
its hopping operators. This relation allows us to analyze 
emerging fermions in three dimensions. We have also 
shown that there is a close connection between strings 
and fermions. The statistical algebra (4) is fundamen- 
tal to this connection since it allows us to determine the 
statistics of the ends of strings from the structure of the 
string operator. 

It is interesting to put this string picture of emerging 
fermions in context. Indeed, it is well-known that gauge 
theories and strings are closely related - the strings corre- 
spond to electric flux lines in the associated gauge theory. 
[23-29] Thus, it appears that the fundamental concepts 
of Fermi statistics, gauge theory, and strings are all con- 
nected. They are all just different aspects of a kind of 
order - topological order - in bosonic lattice systems. 

This research is supported by NSF Grant No. DMR- 
01-23156 and by NSF-MRSEC Grant No. DMR-02- 
13282. 



APPENDIX A: RELATIVE STATISTICS 

In this section, we derive an algebraic formula, analo- 
gous to (4), for the relative statistics of distinguishable 
particles. We begin with a Hilbert space which describes 
two types of hard-core particles hopping on a 2D lat- 
tice. For concreteness, say that there are m particles of 
type 1 and n particles of type 2. The states can be la- 
beled by listing the positions of the two types of particles: 

A typical Hamiltonian for this system is of the form 

<ij> 

where t\j , tj^ are hopping operators for the two types of 
particles. 

We wish to calculate the relative statistics of the two 
types of particles. As before, the statistics are related to 
the simple algebraic properties of the hopping operators. 
Specifically, we will show that the particles have relative 
statistics e'"^ if 

Here, i, k,j, I are (distinct) neighbors of p oriented in the 
clockwise direction. 

We begin with the state which contains a type 

1 particle at z, a type 2 particle at j, and other particles 
which are far away. 

Imagine that we wind particle 1 around particle 2 us- 
ing an appropriate product of hopping operators. One 
expression for the final wound state is: 

wound) = 

{^an---^pq^qi)i^j'a)i^aj)\h j---) 

Here, we've grouped the hopping operators into 6 terms. 
The cumulative effect of these terms is to wind one par- 
ticle around the other (see Fig. 7a) . 
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Similarly, wc can construct a state where the two par- 
ticles don't wind around each other (see Fig. 7b): 

unwound) = 

(4')(4a)(4)K;j-) 

As before, the two states involve the same total path in 
configuration space. Thus, the phase difference between 
the two states is precisely the relative statistics of the 
particles. This can be made rigorous using an argument 
similar to the exchange statistics case. 

We now calculate this phase difference. We use the 
assumed algebraic relation 

Applying this relation to the wound state, and reordering 
the hopping operators using locality, we find 

wound) = e*''^|i;j, unwound) 

This establishes the desired result: the particles have rel- 
ative statistics e*"^. 

A good consistency check for this formula can be ob- 
tained by considering the relative statistics of two parti- 
cles of the same type. In this case, we expect the angle 
for the relative statistics to be exactly twice the angle for 
the exchange statistics: 

g#rel. _ g2iSstat. 

One way to see this is that exchanging the particles twice 
in the same direction is topologically equivalent to wind- 
ing one particle around the other. 

This result, which has a topological character to it, 
can be derived algebraically from our two formulas. We 
start with the expression {tiptpj){ti~ptpi). Applying the 
exchange statistics formula Eq. (4), we find: 

(tiptpj^itkpipl^ — tipitpjtkptpl^ 

— ^zp(^ ipl^kp^pj^ 

Applying the formula again gives 

^ipi^pl^kp^pj) (j^ip^pl^^kp)ipj 

_ tkptpltip^tpj 



Combining these two equations gives 

(tiptpj^itkpipl) —e (tkptpl)(j'ipipj) 

Comparing this with the relative statistics formula 
Eq. (Al), we see that e*"^"!- = e'^^^"*''*- , as claimed. 

APPENDIX B: DIRAC MATRIX 



,ab 



REPRESENTATION OF 7 

In the following, wc represent 7°^'s in terms of the 
Dirac matrices. Note that Eq. (12) implies that 

Thus 7^^, 7^^, 7*'^, and 7^^ satisfy the algebra of Dirac 
matrices. Introducing four Dirac matrices 



7^ =cr^ O a^, 



7^ =a^ (g) 



we can write 



7 =7 , a = x,x,y,y. 



(Bl) 



Similarly, 7^^, 7^^, 7**^, and 7^^ satisfy the algebra of 
Dirac matrices. We can express those operators as 



7°^=i7«75, a = x,x,y,y. 



(B2) 



where 7^ = 'y'^^^'^y^y. Finally, for a,b = x,x,y,y, we 
have 



(B3) 



In this way, we can express all the 7°'', a,b G 
{a;, X, y,y,z,z}, a^h in terms of Dirac matrices. 
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